Chapter 8 

Spin-Orbit Interaction 


8.1 Origin of spin-orbit interaction 

Spin-orbit interaction is a well-known phenomenon that manifests itself in lifting 
the degeneracy of one-electron energy levels in atoms, molecules, and solids. In 
solid-state physics, the nonrelativistic Schrodinger equation is frequently used 
as a first approximation, e.g. in electron band-structure calculations. Without 
relativistic corrections, it leads to doubly-degenerated bands, spin-up and spin- 
down, which can be split by a spin-dependent term in the Hamiltonian. In this 
approach, spin-orbit interaction can be included as a relativistic correction to 
the Schrodinger equation. 

Let us briefly review the origins of this correction, following approach of 
Ref. [80]. For this purpose one has to consider the Dirac equation, which is the 
basic equation for electronic systems, including the electron spin and its rela¬ 
tivistic behavior. One obtains the Dirac equation by linearizing the relativistic 
generalization of the Schrodinger equation. It is Lorentz-invariant and describes 
the electron spin and spin-orbit coupling from first principles. One naturally 
arrives at the Dirac equation when starting from the relativistic expression for 
the kinetic energy 


H 2 = c 2 p 2 + mV. (8.1) 

Inclusion of the electric and magnetic potentials, </> and A, by substituting 
p — (|)A for p (p is the canonical momentum) and H — ecj) for H (electron 
charge e — —e) in equation 8.1 leads to: 

(.H — ecj)) 2 = (cp — sA) 2 + m 2 c 4 . (8.2) 

By interpreting p and H as operators 

d 

p —iKV and H = zft— , (8.3) 

one obtains the relativistic wave equation in external electric and magnetic 
fields. One can write the force-free form of the wave equation as 
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{H 2 -c 2 Y J Pl~m 2 c A )^ = R (8.4) 

P/t = Pxi Pyi Pz are components of the momentum operator, as follows from 
Eq. 8.2. A disadvantage of this differential equation is that ones needs the initial 
values of ^ and to solve it, whereas the Schrodinger equation requires only 
the initial value of The equation can be expressed in the form 

(H - - (3mc 2 )(H + c^a^p^ - (3mc 2 )ip = 0 (8.5) 

A' /-' 

if the constant coefficients and /3 satisfy the relations 

a^f3 + f3a li = 0 
/ 3 2 = 1 . 

One can linearize it by taking only the first part of Eq. 8.5. Solution of the 
first part will give the solution of whole equation. The linearized equation has 
the advantage that it is of first order in similar to the Schrodinger equation; 
it is referred to as Dirac equation: 

(H - cY^otyPn - f3mc 2 )-ip = 0, (8.6) 

with the four-component vector ip. 

Let us now compare the Dirac equation with the Schrodinger equation. For 
this we use the non- linearized form 


[H — £(/) — ca • (p—-A) — /3mc 2 ][H — £(f) + ca • ( p — - A) + = 0. (8.7) 

c c 

Using the approximation that the kinetic and potential energies are small com¬ 
pared to me 2 , two components of the spin function can be neglected, and equa¬ 
tion 8.7 takes the form 


^- £ -c.N 
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Am 2 c 2 


cr*(E x p) ] f = W'tjj, 

(8.8) 


where W + me 2 is the total energy. 

The first two terms in the parentheses are equivalent to those in the Schrodinger 
equation for external fields. The third term corresponds to the interaction en¬ 
ergy —jl • B of the magnetic dipole, whose moment is presented by operator 
/2 = 2 ^^ = ^s. The fourth term is a relativistic correction to the energy 
and does not have a classical analogy. The fifth term describes the spin-orbit 
coupling. 

It can be illustrated in framework of the classical electrodynamics, where 
the vectors of the electromagnetic field depend on the reference system. In a 
reference system that moves with velocity v relative to an electric field E, one 
finds a magnetic field 
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B =--v x E = —(E x p), (8.9) 

c me 

where terms of order (^) 2 and higher order terms are neglected. In other words, 
the moving electron experiences a magnetic field in its rest frame that arises 
from the Lorentz transformation of the static (external) electric field; this field 
will affect the electron spin. The energy of the electron in this field , due to its 
magnetic moment /i, is 

— fi • B = ——— s • B =-^-~s-(E x p). (8.10) 

me m z c z 

This additional energy term in the Hamiltonian is essentially the spin-orbit 
term in Eq. 8.8, except for a factor of 2. This factor (Thomas factor) is missing 
because we have not taken into account that changing the frame of reference 
also leads to a time transformation and, consequently, the precession frequency 
of the electron spin in the magnetic field changes. 

In case of centrally symmetric electrical fields, as e.g. the orbital motion of 
an electron in the electric field of an atomic nucleus, one has 


E = 


lr dV_ 

er dr ’ 


and the term can be written in the form (s = |cr) 


eh 


Am 2 c 2 


cr-(E X p) = 


_lrdK 
2m 2 c 2 S l e r dr 
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P = 


2 m 2 c 2 r dr 


IdV 


( 8 . 11 ) 


( 8 . 12 ) 


This is called spin-orbit energy, as it results from an interaction of the spin 
with the magnetic field that is experienced by the moving electron. 


8.2 Inversion symmetry and time reversal symmetry 

The form of the spin-orbit interaction term in Eq. 8.2 suggests that a spin- 
degenerate level can be split by spin-orbit coupling (SOC) into levels with spin 
parallel and antiparallel to the orbit. In a free atom, spin-orbit interaction can 
lift the degeneracy of states with the same orbital (spatial) wave function but 
with opposite spins. In solids, however, such a splitting can be forbidden due 
to crystal symmetry. 

As a fundamental principle of quantum mechanics, time reversal symmetry 
always preserves the Kramer’s degeneracy between a function ^(r, s) and its 
complex conjugate / 0*(r, s); the latter function differs from ^ through a simul¬ 
taneous reversal of wave vector and electron spin. Therefore, at any point of 
the Brillouin zone, one can write for the energy 

E(k,r)=E(-k,i). (8.13) 

The Kramer’s degeneracy remains unaffected when including the spin-orbit 
interaction term in the Hamiltonian. (Notice that j and j do not necessar¬ 
ily denote eigenfunctions of the spin operator a z , but correspond to the two 
quantum numbers denoting eigenstates for a given Bloch vector.) 
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If the crystal lattice has inversion symmetry (i.e. if the operation r —> —r 
does not change the crystal lattice), one will obtain 

£(M) =£(-M) and E(k,l)=E(-k, |). (8.14) 

From combination of the equations 8.13 and 8.14 it becomes clear that if both 
time reversal symmetry and inversion symmetry are present, the band structure 
should satisfy to the condition 


E(k,]) = (8.15) 

In other words, the energy cannot depend on the electron spin. Conse¬ 
quently, for crystals which have inversion symmetry (like fee, hep), spin splitting 
is not allowed in the bulk, and these solids keep their spin degeneracy. 


8.3 Hydrogen atom 

Consideration of the electron orbital motion in the Coulomb field permits us 
to visualize spin-orbit interaction. We use the hydrogen atom for which the 
Schrodinger equation can be solved analytically. In the semiclassical model, 
moments are treated as vectors with z-components being already quantized. 
(For spin-orbit interaction, only the case of l > 1 is important, since SOC 
vanishes for 1 = 0.) The energy-level correction due to spin-orbit coupling takes 
the well-known form [81]: 

En,i,s — E n — fi s • Bi = E n + g^° m 2 r 3 ( s ' 1 )’ (8.16) 

where 


fio Ze ^Ze 

Bi = ——(v x r) = 1 . 


(8.17) 


47r 3 47r 3 m 

The dot product of spin and orbital momentum can be expressed though 
the eigenvalues of spin, orbital, and total angular momentum 


s • 1 — -[j 2 - l 2 - s 2 ] — -h 2 [j(j + 1 ) - l{l + 1 ) - s(s + 1 )]. ( 8 . 18 ) 

Therefore the energy can be written as 

E n ,i,s — E n + -\j(j + 1 ) — 1(1 + 1 ) — 5(5 + 1 )], ( 8 . 19 ) 

where a = §^ 2^3 ? - the so-called spin-orbit coupling constant. Thus, the 
energy level defined by n, l splits into sublevels with different total moments, 
j = l + \ and j = l — 

In quantum mechanics, the electron radius is treated using the probability 
density distribution. So the observable of \ should be written as its expectation 
value 
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a = 


/ ioZe 2 h 2 




StCTU 2 / ^n,l,m r3 ^n,l,mdT. 


( 8 . 20 ) 


If we substitute the Hydrogen wave function for / we obtain 


CL — E n 


Z 2 a 2 

n • 1(1 + 1/2) (Z + 1) ! 


( 8 . 21 ) 


2 2 i 

with the relativistic fine-structure constant a = — jgf- Finally, 

the energy difference between two sublevels with opposite spin is 


Z 2 a 2 


AE/ S = —E n 


n • /(/ ± 1) 


( 8 . 22 ) 


From Eq. 8.22 we immediately see two important aspects of SOC: (1) atomic 
spin-orbit splitting is proportional to the product of E n -Z 2 ; in case of a Coulomb 
field, SOC is proportional to (2) Spin-orbit splitting decreases with increas¬ 
ing quantum numbers n and l . 


8.4 Solids 

Bulk crystals 

If a crystal lattice potential has a center of inversion, spin-orbit interaction will 
not affect the band structure, except for certain points of the Brillouin zone 
(BZ), where the influence of spin-orbit interaction can be very important. As 
an example, we consider the center of a BZ (r point) that has cubic symmetry. 
In a tight binding model, we can build up the p bands from atomic p wave 
functions . In a free atom without SOC, there are three degenerate p functions, 
so three p bands will appear in the model . These bands will all be degenerate 
at the T point (k — 0, see Fig. 8.1(a)), transforming one into the other through 
cubic symmetry transformations at this point. Each of the three bands is also 
doubly degenerate in the spin, altogether a six-fold degeneracy at T. 

When spin-orbit interaction is included, the bands will split in two sub¬ 
bands with fourfold p 3 / 2 and twofold p 1 / 2 degeneracy (see Fig. 8.1(b)). The 
bands will now behave differently as we move away from the V point depending 
on whether or not there is inversion symmetry in the lattice. If the crystal has 
a center of inversion, each band will preserve its spin degeneracy, as illustrated 
on Fig. 8.1(b). The p 3 / 2 bands, with j — 3/2, will now split into two bands 
with rrij = ±3/2 and rrij = ±1/2. (with projections taken e.g. in the direction 
of the vector k). 

If a center of inversion symmetry is absent, then the p 3 / 2 and pi/ 2 bands will 
additionally be split by spin-orbit interaction, removing the spin degeneracy. 
This situation is shown in Fig. 8.1(c). 

The cases considered are really occurring in crystals, e.g. in semiconductors 
like Ge and InSb. Germanium has a center of inversion symmetry, and we will 
have the first scenario for the top valence band. For InSb, which crystallizes in 
a zinc-blende lattice, the spin degeneracy will be lifted as we leave the center of 
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Figure 8.1: Influence of spin-orbit interaction on the p-levels at the center of 
the Brillouin zone: (a) Without SOC, six degenerate levels are observed at the 
T-point; (b) SOC leads to a splitting into p 3 / 2 and pi/ 2 levels, but leaves the 
spin degeneracy in lattices with inversion symmetry; (c) in a lattice without 
inversion symmetry, the spin degeneracy is fully removed, except for the T- 
point. 


the BZ. This effect of symmetry on the electronic band structure of zinc-blende 
crystals leads to the appearance of a spin-orbit interaction term that is cubic 
in k (Dresselhaus splitting). 

Surfaces 

Termination of the crystal by a surface breaks the 3D inversion symmetry. In 
this way, it permits electron surface states with the same parallel wave vector, 
but fey, opposite spins to have different energies: 

E(k h ])^E(k h i). (8.23) 

In the absence of external or internal magnetic fields, spin-orbit splitting is 
the only possible cause for a splitting of spin degenerate levels at the surface, 
and this interaction is described by the Rashba Hamiltonian [82]. 













